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A  BSTRACT 


The  acoustic  impedance  of  several  Helmholtz  resonators 
were  measured  using  two  different  techniques.  The  first 
technique  employed  the  conventional  Standing  Wave  Tube 
method.  The  second  was  the  Dual-Microphone  Transfer  Function 
method  of  J.T.  Chung  and  D.  A.  fllaser.  The  calculability  of 
the  frequency  dependent  Helmholtz  resonator  impedance 
allowed  the  methods  to  be  compared  for  both  absolute  accu¬ 
racy  and  relative  precision.  The  average  relative  precision 
for  the  SWT  technique  was  4.2X  while  the  average  for  the 
dual-microphone  technique  was  18.4SS. 
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i.  introduction 

acoustic  impedance  tubes  have  been  used  for  many  years 
to  obtain  complex  acoustic  impedances,  reflection  coeffi¬ 
cients  and  phase  information  for  nar.y  different  types  of 
materials.  The  techniques  used  with  the  impedance  tubes  have 
been  many  and  varied.  This  study  involves  the  comparison  of 
the  dual-microphone  transfer  function  technique  used  by 
Chung  and  Blaser  [Ref.  1]  and  the  standing  wave  tube  (SWT) 
technique  to  determine  their  absolute  accuracy  and  relative 
precision.  This  work  was  undertaken  as  a  prelude  to  the 
extension  of  computer  controlled  impedance  measurement  tech¬ 
niques  to  water. 

A.  IHPEDANCE  MEASUREMENT  TECHNIQUES 

Many  different  impedance  measurement  techniques  have 
been  used  over  the  years.  Baranek  [Ref.  2]  divides  the 
different  techniques  into  three  general  categories.  Surface 
Methods,  Transmission  Line  Methods  and  Comparison  Methods. 
This  study  looked  closely  at  only  two  of  the  techniques, 
however,  seme  of  the  other  methods  are  worth  mentioning  as 
background. 

1  •  Surface  Methods 

These  methods  involve  the  measurement  of  the  sound 
pressure  at  a  point  on  the  surface  of  interest  and  the 
particle  or  volume  velocity  at  the  surface.  Cook  [Ref.  3] 
did  this  using  a  short  tube  of  a  given  radius  and  length  and 
a  clamped-edge  diaphram.  He  was  then  able  to  calculate  the 
volume  velocity  and  the  average  pressure  over  the  surface  of 
the  sample.  Beranek  [Ref.  « ]  presented  a  useful  theory  for 
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predicting  the  normal  acoustic  impedance  of  porous  material? 
in  terms  cf  their  flow  resistance,  porosity  ar.a  err ec*:v- 
dynaaic  mass  of  the  enclosed  gas. 

2.  Transmission  Line  M  ethods 

By  far,  the  most  popular  methods  of  impedance  meas¬ 
urements  fall  into  this  category.  The  standing  wave  tube  and 
dual-microphone  techniques  in  this  study  are  included.  They 
will  be  examined  in  detail  later.  Another  technique  involves 
the  measurement  of  the  pressure-frequency  resonance  curve 
and  determining  the  curve  width  and  resonant  frequency. 
Harris  [Ref.  5]  [Ref.  6]  discusses  two  of  these  methods 
which  are  closely  related  to  each  other. 

Finally,  one  additional  method  currently  in  use  by 
the  O.S.Navy  for  making  acoustical  measurements  deserves 
mentioning.  [Ref.  7]  This  system  is  digital  and  is  used  for 
underwater  acoustical  measurements  from  1  Hz.  to  approxi¬ 
mately  2  M  Hz.  This  digital  measuring  system  employes  a 
buffered  video  analog-to-di gital  converter  and  was  primarily 
designed  to  perform  acoustical  measurements  using  the  Prcny 
Method.  [Ref.  8]  [Ref.  9]  The  system  uses  a  stable  signal 
and  sampling  source  and  12-bit,  5  MHz.,  analog-tc-digital 

converter  with  a  4K  buffer  memory.  A  PDP-11/23  serves  as  a 
controller  as  well  as  a  signal  processor  for  averaging  and 
amplitude  and  phase  estimation. 

3 •  Comparison  Methods 

These  methods  mainly  involve  the  use  of  acoustical 
bridges  or  reliable  acoustical  standards.  Steward  first 
introduced  [Ref.  10]  the  simplest  and  most  rapid  to  use. 
Later  refinements  followed  [Ref.  11]  but  suitable  variable 
standards  of  acoustical  resistance  and  reactance  were  diffi¬ 
cult  to  obtain.  The  use  cf  reliable  acoustic  standards, 
rigid  wall  and  an  eighth-wa velength  tube,  was  reported  by 
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Flanders  (Ref.  12]  in  1932.  HoksV®;,  mathematical  calcula¬ 
tions  were  difficult. 

In  this  study,  the  transmission  line  method  an* 
comparison  method  are  combined  by  introducing  a  r.ev  "  stan¬ 
dard"  #  the  Helmholtz  resonator.  This  standard  has  the 
advantage  of  having  an  easily  calculable  reactance  which  is 
absolute  to  within  one  adjustable  parameter.  This  parameter 
is  fixed  by  the  Helmholtz  resonance  frequency  which  is 
easily  measured  to  a  high  degree  of  precision  and  absolute 
accuracy . 

B.  HISTORY  OF  DUAL- MICROPHONE  MEASURING  TECHNIQUES 

ClapF  and  Firestone  [Ref.  13]  first  mentioned  the  use  of 
two  microphones  in  an  article  published  in  19# 1.  In  this 
article  it  was  recognized  tha*  their  equipment  could  be  used 
tc  measure  the  absorption  coefficient  of  materials.  Later, 
Shultz  [Hef.  14]  published  an  article  on  the  acoustic  watt¬ 
meter  and  its  application  for'  measuring  absorption  coeffi¬ 
cients.  During  this  time  the  use  of  the  acoustic  wattmeter 
fcr  impedance  measurements  was  little  used  since  the 
standing  wave  tuba  (SWT)  technique  was  the  most  favcred 
method.. 

In  1977  Seybert  and  Ross  [Ref.  15]  published  an  article 
which  discussed  an  experimental  technique  using  two  micro¬ 
phones  and  random  noise  excitation.  This  technique  used 
digital  methods  and  the  aut  c-  and  cross- spectra  to  determine 
the  normal  complex  reflection  coefficient  cf  a  sample  termi¬ 
nating  a  tube.  Finally,  in  1980,  Ihung  and  Blaser  [Ref.  1] 
published  their  article  discussing  the  transfer  function 
technique  fcr  measuring  in-duct  acoustic  properties.  This 
technique  also  used  digital  methods  and  the  auto-  and 
cross-spectra  to  determine  the  normal  complex  reflection 
coefficient.  However,  the  mathematics  did  r.ct  involve  the 

1 1 


matrix-inversion  used 


Ross  and  as 


by  Seybert  aid 
the  computation  time  was  significantly  reduced. 

A  two  microphone  technique  utilising  pure 
than  random  excitation  ari  analog  rather 
processing  is  described  by  Elliot.  [Ref.  16] 


C. 


C.  POHPCSE  OF  STUDY 

The  above  methods  have  been  used  to  measure  numerous 
materials  of  differing  acoustic  properties.  However,  the 
properties  of  these  materials  were  not  known  or  calculable 
to  any  reasonable  accuracy.  The  aim  of  this  study  was  to 
determine  the  absolute  accuracy  of  the  dual-micrcphcne  tech¬ 
nique  and  its  relative  precision  compared  to  the  SWT  tech¬ 
nique  as  a  prelude  to  the  development  cf  an  impedance 
measurement  apparatus  for  use  in  water.  A  Helkholtz  reso¬ 
nator  was  used  to  accomplish  this  comparison. 
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II.  COMPARISON  OF  MEASUREMENT  TECHNIQUES 

Previous  studies  of  impedance  measurement  techniques 
involved  the  use  of  materials  for  which  the  impedance  was 
net  accurately  known.  This  study  deals  with  impedance  meas¬ 
urements  using  a  Helmholtz  resonator  as  a  calculable  stan¬ 
dard  to  determine  the  absolute  accuracy  and  relative 
precision  of  the  two  techniques.  A  Helmholtz  resonator  was 
chosen  because  its  various  acoustical  parameters  could  be 
independently  measured  or  calculated  and  then  compared  to 
the  experimental  results. 

A.  THEORETICAL  BACKGROUND 

The  term  "acoustic  impedance"  can  have  different  mean-, 
ir.gs  for  different  situations.  It  is  necessary  to  affine 
exactly  the  different  impedances  involved  in  this  study. 
This  is  necessary  to  prevent  confusion  later  when  the  rela¬ 
tionship  between  the  different  impedances  is  used  to  deter¬ 
mine  the  calculated  acoustical  properties  of  a  Helmholtz 
resonator.  Normally  this  is  not  necessary  since  the  specific 
acoustic  impedance  and  the  acoustic  impedance  differ  by  a 
constant  with  the  dimensions  of  area. 

1  •  Specific  Aco  ustic  I  mpe dance 

The  specific  acoustic  impedance  is  defined  as  the 
ratio  of  acoustic  pressure  in  a  medium  to  the  associated 
Eulerian  particle  speed.  [Ref.  17] 

Zs  *  p/u  (2.1) 

The  MKS  unit  of  specific  acoustic  impedance  is  ?a-s/m. 
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The  specific  acoustic  impedance  of  a  medium  is  a 
real  quantity  for  plane  waves.  However,  this  net  tru-  for 
the  standing  waves  generated  in  an  impedance  tube.  For  this 
case,  the  specific  acoustic  impedance  will  be  complex  and 
can  expressed  in  the  form 

Zj  =  r  +  jx  (2.2) 

where  "r"  is  the  specific  acoustic  resistance  and  "x"  is  the 
specific  acoustic  reactance  for  the  medium  concerned. 

2.  Acoustic  Impedance 

The  acoustic  impedance  of  a  medium  acting  cn  a 
surface  cf  area  "A”  is  the  ratio  of  the  acoustic  pressure  at 
the  surface  to  the  volume  velocity  at  the  surface  [Ref.  17] 

Zo.  =  p/0  (2.3) 

where  0  =  uA  if  u  is  normal  to  the  surface  of  area  "A”.  The 
MKS  unit  of  acoustic  impedance  is  Pa-s/m^,  which  is  also 
termed  an  acoustic  chm.  The  acoustic  impedance  and  the 
specific  acoustic  impedance  are  related  at  the  surface  of 
area  "A"  by 

Zx  =  Zj  /  A  (2.4) 


3.  Mechanical  Impedance 

The  mechanical  impedance  is  the  ratio  cf  the  driving 
force  to  the  resultant  particle  speed  at  the  point  where  the 
force  is  applied  [Ref.  17] 

Z*,  =  F/u  (2.5) 

The  HKS  unit  of  mechanical  impedance  is  N-s/m,  also  referred 
to  as  a  mechanical  ohm. 
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Helmholtz  5 3 senator 

A  Helmholtz  resonator  can  be  treated  as  a  lumped 
parameter  system.  The  resonator  consists  of  a  volume  "7" 
with  rigid  vails  and  a  neck  of  length  "L"  and  cross- 
sectional  area  "A".  The  acoustic  pressure  in  the  volume 
provides  the  lumped  stiffness  "s'*  if  the  wavelength  is  much 
greater  than  the  cube  root  of  the  volume.  The  fluid  in  the 
neck  provides  the  lumped  mass  as  long  as  the  wavelength  is 
much  greater  than  the  length  of  the  neck.  Finally,  the 
resistance  is  provided  by  the  opening  radiating  sound,  if 
the  wavelength  is  much  greater  than  the  square  root  of  the 
cross-sectional  area  of  the  neck.  An  additional  resistance 
is  introduced  by  visccus  losses  in  the  neck  and  irreversible 
thermal  conduction  to  the  volume  boundries  if  the  compres¬ 
sion  of  the  gas  is  adiabatic  and  the  walls  remain  isoth¬ 
ermal.  However,  since  these  resistance  sources  are  not 
easily  calculated  from  '‘First  Principles",  they  are 
neglected.  The  radiation  resistance  "R  "  can  be  calculated, 
however,  its  value  is  dependent  upon  the  impedance  of  the 
medium  which  changes.  To  get  the  viscous  loss  resistance 
term  the  flew  field  in  the  neck  is  required.  Therefore  these 
resistance  terms  are  of  little  use  in  this  calculation. 

The  fluid  in  the  neck  has  an  effective  mass  given  by 
[Ref.  17] 

m  =  /o  AL '  (2.6) 

where  "L'  "  is  the  effective  length  of  the  neck  and  is  a 

function  of  actual  length  "L"  and  the  neck  radius  "a"  and 
the  flew  field  in  the  vicinity  of  the  oriface.  The  effective 
mass  is  greater  than  the  actual  mass  in  the  neck  due  the 
radiation-mass  loading  at  the  ends  generated  by  this 
extended  flow  field.  The  stiffness  of  the  resonator  is 
given  by  [Ref.  17] 

s  *  (yOC*)  (A*  )/7  (2.7) 
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By  driving  the  Helmholtz  resonator  with  a  sound  wave 
of  pressure  amplitude  "P",  the  following  equation  for  Th¬ 
in  ward  displacement  "x"  of  the  fluid  in  the  r.eck  results 
[Ref.  17] 

nyjk  +  *  sx  ~  APeJ**  (2.8) 


This  differental  equation  is  the  same  as  that  for  a  driven 
oscillator  and  the  mechanical  impedance  of  the  Helmholtz 
resonator  can  be  expressed  as  [Ref.  17] 

Z„  *  Rr  ♦  j  (*>m  -  s/j  )  (2.9) 

where  the  mechanical  reactance  '‘X,*"  is  given  by 

ZA  -  sfo  (2.  10) 

Since  the  resonance  occurs  when  X*,  *  0,  the  resonant 

frequency  "  follows 

uj,1  *  s/m  (2.  11) 

Using  equations  (2.6)  and  (2.7) 


The  final  result  obtained  from  cancelling  terms  is 


tOa  - 


(2. 12) 


It  should  be  noted  that  nowhere  in  the  derivarion  of 
equation  (2.12)  is  any  restriction  placed  on  the  physical 
arrangement  of  the  Helmholtz  resonator.  Figure  2.1  shows  the 
arrangement  used  for  this  study.  As  can  be  seen  from 


16 


equation  (2.12)  and  the  figure,  the  resonant  frequency  car. 
ta  chanced  by  changing  the  area  of  «!ieckM  (diameter  T>,  \  , 
volume  of  the  resonator  (change  depth  h)  or  the  length  of 
!?,  n  r*  ^  *  r  i.  c  ,vc  r.®  ^  s  * ) 

Recalling  equation  (2.10) 

X„  =  tdm  -  s4j  (2.  10) 

and  fcy  factoring  out  "a"  and  using  equation  (2.11),  the 
following  result  is  obtained 

xm  s  m  ( ) 

Since  u>  =  2>7*f  the  final  form  is  obtained 
X„  =  2m»(f  -  f.v/f) 

where  f  is  the  frequency  in  hertz.  For  convenience  we  define 
X*  »  f  -  f|/f  (2.  13) 

so  that 

X„  *  23rmXi,  (2.14) 

It  is  now  possible  to  derive  the  relationship 
between  the  specific  acoustic  reactance  "X,  '•  and  "Xi".  Given 
a  resonator  as  shown  in  Figure  2.2,  divided  into  sections  1 
and  2  as  indicated  and  recalling  equation  (2.1) 

Zs  =  p,  /u,  (2.15) 

where  d/  and  u,  represent  the  pressure  and  particle  speed, 
respectively,  in  section  1.  This  is  the  specific  acoustic 
impedance  of  the  medium  in  section  1. 

Next  using  equation  (2.5) 

Z„  *  (2.16) 

Converting  F*  zo  a  pressure  at  the  boundary  results  in 


I 


Z, 


(2. 17) 


•»  ut 


where  a,  and  u,  represent  the  Dressure  and  partrcis  steed, 
respectively,  in  section  2.  A*  is  the  cross-sectic r.al  area 
of  the  resonator  neck.  Since  there  is  no  physical  boundary 
between  sections  1  and  2,  the  pressures  on  both  sides  must 
be  equal. 

P,  3  Fz  *  P 

Also,  the  volume  velocities  must  be  equal  across  the 
boundry.  Therefore 

A,  u,  =  AJtui 

ut  =  yjk-u,  (2.18) 


and 


Substituting  equation 
regrouping  the  appropriate 


(2.18)  into 
terms  results 


equation 

in 


(2.  17) 


z*>  = 


(2.  19) 


Finally,  it  can  be  seen  that  the  quantity  (p/u,)  is  just  the 
specific  acoustic  impedance  and  the  result  is 


(2.20) 


Since  the  impedances  are  composed  of  the  resistive  and  reac¬ 
tive  parts,  the  reactive  portions  vary  in  the  same  manner. 
Therefore, 


(2.21) 
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It  can  be  seen  from  equation  (2 . 1  a >  that 


=  2  ZTm-jfc- 


It  should  be  noted  that  the  two  techniques  in  this 
study  provide  the  values  for  specific  acoustic  impedance 
normalized  to  the  characteristic  impedance  n />  c" .  Therefore, 
the  preceding  equation  becomes 


X*,  A*, 


and,  finally,  substituting  for  ”m'*  from  equation  (2.6)  and 
for  "X5/X*,  "  from  equation  (2.21) 


_ 


(2.22) 


As  can  be  seen,  this  equation  does  not  depend  on 
frequency  and  therefore  remains  constant  as  frequency 
changes.  For  convenience  this  constant  is  called  "K*/'  with 
the  subscript  signifying  the  theoretical  value  for  the  ratio 
of  "X, /fic'C,"  for  a  given  Helmholtz  resonator. 


jlza'A _ 

C/4. 


(2.23) 


"  has  been  reduced  to  an  expression  which 
contains  only  directly  measurable  terms  with  the  exception 
of  ” L ,  the  effective  length  of  the  neck.  However,  by  using 
equation  (2.12)  ML,M  can  be  determined  from  parameters  which 
are  measurable  by  separate  means.  Later  in  this  chapter, 
"L* "  will  be  looked  at  more  closely  for  comparison  with  the 
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calculated  "L"*  for  ether  Helmholtz  resonator 
conf i gurat ions . 

B.  STANDING  WAVE  TUEE  TECHNIQUE  (SWT) 

The  standing  wave  tube  technique  provides  a  convenient 
laboratory  method  for  determining  the  absorption  coefficient 
ar.d  specific  acoustic  impedance  of  materials.  The  method, 
however,  is  somewhat  tedious  when  a  large  number  of  samples 
is  measured  or  measurements  of  a  single  sample  are  required 
at  a  great  number  of  frequencies. 

1 .  Theory 

The  sample  to  be  measured  is  placed  at  one  end  of 
the  standing  wave  tube.  At  the  other  is  mounted  a  loud¬ 
speaker.  A  microphone  probe  is  inserted  through  the  loud¬ 
speaker  end  and  is  attached  to  a  microphone  car.  The  car 
moves  such  that  the  microphone  probe  is  moved  along  the  axis 
of  the  tube. 

A  pure  tone  signal  is  reflected  from  the  sample.  The 
reflected  wave  interferes  with  the  incident  wave  and  a 
standing  wave  is  established.  The  microphone  probe  and  car 
arrangement  allows  fer  the  location  and  acoustic  pressure  of 
the  nodes  and  anti-nodes  to  b9  measured.  From  the  measure¬ 
ment  of  these  amplitude  and  phase  characteristics,  calcula¬ 
tion  of  absorption  coefficient  and  specific  acoustic 
impedance  of  a  sample  is  possible. 

From  reference  6,  the  following  are  the  equations 
needed  for  the  determination  of  the  desired  acoustic 
parameters. 


X 

>~c 


IfU. . 


where 


(2.24) 
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(2.25) 


> 


^  =  cot  h 


'[-fef-] 


and 

(2.26) 

and  minimum  pressures  in  the 
the  distance  to  the  first 
pressure  minimum  (node)  from  the  sample.  (NOTE:  The  micro¬ 
phone  output  voltage  is  the  actual  measurement  made  since 
only  the  ratio  of  the  pressures  is  required) 

2 .  Equipment  Set-up 

Figure  2. 3  is  a  simplified  representation  of  the 

equipment  used  in  the  standing  wave  tube  technique.  It 

includes  the  Bruel  and  Kjaer  Standing  Nave  Apparatus  Type 
400  2,  a  HP  302A  Wave  Analyzer  and  an  amplifier  and  frequency 
counter. 

a.  Bruel  and  Kjaer  Type  4002  Standing  Wave 

Appa  ratus 

The  apparatus  is  designed  for  easy  and  quick 
evaluation  of  acoustic  materials  by  the  standing  wave 
method.  A  major  advantage  is  that  relatively  small  circular 
samples  (10  cm.  in  diameter)  are  all  that  are  required. 

The  tube  used  had  an  interior  diameter  of  10 
cm.,  which  covered  a  frequency  range  of  90  to  1800  kz.  The 

lower  limit  was  set  by  the  length  of  the  tube  and  the  upper 

limit  determined  by  the  cut-off  frequency  for  the  lowest 
non-plane  wavemode  of  the  tube.  [Ref.  18]  Two  sample 

holders  were  also  used  with  fixed  depths  of  1  in.  and  2  in., 
respectively.  The  bases  of  the  holders  were  very  thick  to 
present  as  rigid  a  termination  as  possible. 
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The  microphone  in  the  microphone  car  is  a  snail 
dynamic  microphone  of  adequate  frequency  range  for 
standing  wave  tube  method.  It  is  isolated  by  elastic 
supports  in  the  car  to  reduce  the  effects  of  external 
airborne  noise  and  impact  sound  or  vibration. 

b.  Hewlett-Packard  HP  302A  Wave  Analyzer 

The  Hewlett-Packard  model  302A  wave  analyzer  is 
a  tunable  voltmeter  cf  high  sensitivity  over  the  frequency 
range  of  20  Hz.  to  5 0  kHz.  The  wave  analyzer  in  the  BFO  mods 
of  operation  becomes  an  audio  signal  generator  combined  with 
a  tuned  voltmeter.  Further,  since  the  tuned  voltmeter  has  a 
very  narrow  pass  band  (6  Hz) ,  noise  and  harmonics  of  the 
signal  have  very  little  effect  on  the  voltmeter  reading. 

3 .  Experimental  Proced ures 

The  standing  wave  tube  technique  was  used  in  several 
different  experiments.  It  was  used  with  a  rigid  end  as  th® 
sample  to  measure  the  "acoustic"  end  of  the  microphone  probe 
tube  and  the  sound  speed.  It  was  also  used  to  determine  the 
resonant  frequency  and  reactive  impedance  of  three  Helmholtz 
resonators.  Taole  I  shows  the  physical  dimensions  of  each  of 
the  resonators. 

a.  Data  Collection 

In  all  of  the  experiments  the  data  was  collected 
in  the  same  manner.  However,  in  the  first  experiment  (rigid 
end)  the  only  data  collected  was  the  location  of  the  first 
two  nodes  and  the  frequency  of  the  signal.  Table  II  shows 
this  data. 

For  the  experiments  involving  the  Helmholtz 
resonators  the  values  of  the  maximum  and  minimum  microphone 
output  voltages  were  obtained  in  addition  to  the  frequency 
and  node  locations.  This  data  was  placed  in  Tables  III-V 
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T  ABLE  I  i 

I 

Helmholtz  Fesonator  Physical  Dimensions  I 


Fescnator 

Vo  la  ue 

N9£jc  Diameter 

Keck  Length 

no. 

cm . 

cm . 

cm . 

1 

195. 52 

2.52 

1.27 

2 

381 . 01 

2.52 

1.27 

3 

195. 52 

3.46 

1.00 

b.  Data  Seduction 


In  the  first  experiment,  the  measured  distance 
between  the  two  nodes  was  used  to  calculate  the  wavelength. 
From  this,  the  quarter- wavelength  was  determined  and 
compared  to  the  distance  from  the  sample  to  the  first  node 
as  measured  by  the  microphone  car.  Since  the  "sample”  was  a 
rigid  end  located  5.00  cm.  from  the  measured  end  of  the 
tube,  the  two  distances  should  have  differed  by  5  cm. 
Additionally,  the  wavelength  and  the  frequency  were  used  to 
calculate  the  sound  speed  in  the  tube  for  each  frequency 
measured.  These  results  are  shown  in  Table  II  The  average 
value  of  M  d"  suggests  that  the  "acoustic  center"  of  the 
probe  tube  is  1.7  +  .5  mm.  ahead  of  the  physical  end  of  the 
probe  tube  over  the  frequency  range  examined. 

For  the  experiments  involving  the  Helmholtz 
resonators,  the  reduction  of  the  data  was  somewhat 
different.  The  wavelength  was  determined  by  one  of  twc 
methods  depending  upon  whether  two  node  locations  were 
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T  ABLE  II 

Standing  Wave  Tube  Sound  Speed  Data 


! 


Frea . 

First  Min. 

Sec. 

Hz. 

cm. 

cm 

299 

23.90 

81. 

400 

16.80 

59. 

600 

9.58 

38. 

800 

5.96 

27. 

1001 

3  .84 

21. 

1  197 

2.19 

16. 

1402 

1  .40 

13. 

1599 

0.54 

11. 

Min 

ad 

c 

• 

cm . 

cm . 

cm.  /s 

74 

28.93 

5.03 

34588 

95 

21.53 

4.78 

34520 

38 

14. 40 

4.82 

34560 

53 

10.78 

4.79 

345  12 

03 

8.60 

4.76 

344  14 

91 

7.35 

5.  16 

35240 

64 

6  .  13 

4.73 

34321 

36 

5.40 

4.86 

34602 

average 

4.  87 

34595 

deviation 

.150 

277 

per  cent 

03.08 

0. 80 

J 


measurable.  In  those  cases  where  two  or  more  node  locations 
were  measurable,  the  distance  between  them  was  used  to 
calculate  the  wavelength.  At  the  lower  frequencies  where 
only  one  node  was  measurable,  the  wavelength  was  determined 
from  the  frequency  and  the  average  sound  speed  found  during 
the  experiment.  It  should  be  noted  that  an  attempt  was  made 
to  use  the  location  cf  the  measurable  anti-node  to  determine 
wavelength.  However,  this  measurement  was  too  inaccurate  due 
to  the  nature  of  the  waveform  being  measured.  At  the 
maximum,  the  values  of  amplitude  are  changing  very  slowly 
(second  order)  so  that  the  measurement  of  the  exact  location 
of  the  maximum  was  very  difficult.  On  the  other  hand,  at  the 
nodes,  the  values  were  changing  more  rapidly  {first  order) 
and  the  location  could  be  easily  and  precisely  determined. 


|  TABLE  III 

t  ' 

(  Standing  Wave  Tube  Measured  Data  -  Resonator  #1  I 


Free. 

Fir  s  t 

Min . 

Sec . 

HZ. 

cm . 

uv . 

cm . 

250 

31.75 

65 

350 

20.  35 

50 

69.  70 

400 

16.6  5 

26  0 

59.  69 

450 

12.30 

14  0 

50.  40 

500 

6.  70 

450 

41.  20 

550 

1.70 

240 

30.  35 

600 

.55 

200 

22.  50 

650 

.  10 

200 

18. 42 

700 

15.  90 

750 

14.  20 

1  000 

9.75 

1250 

7.  60 

1  500 

6. 75 

Min. 

Third 

Min. 

Maximum 

uv. 

cm . 

uv. 

mv . 

2.20 

60 

1.30 

260 

4.  80 

130 

39.20 

175 

1.80 

460 

75.30 

520 

2.80 

290 

61.65 

190 

1.60 

180 

51 .  30 

185 

1.78 

140 

45.  10 

155 

2.60 

34 

40. 64 

31 

1.00 

90 

37.  15 

1  00 

5.20 

30 

27.  15 

20 

3.40 

5 

21.40 

5 

1.95 

25 

17.  80 

5 

2.50 

TABLE  IT 

Standing  Wave  Tube  Measured  Data  -  Resonator  #2 


Freq. 

lirsi. 

Min . 

Sec . 

£i£- 

Third 

Min. 

Maximum 

HZ. 

cm . 

uv . 

cm  . 

uv. 

cm . 

u  V. 

mv. 

250 

27.88 

1750 

22.  1 

300 

19.69 

335 

76.  72 

330 

2.90 

350 

10.12 

1610 

59.  32 

1  870 

8.45 

375 

3.65 

1280 

49.  63 

15  10 

6.70 

400 

1.59 

800 

40.  94 

915 

84 .65 

355 

4.35 

4  50 

.43 

480 

29.  90 

615 

68.58 

5  65 

8.00 

500 

0.00 

164 

23.  96 

75 

58.64 

128 

2.88 

750 

13.53 

24 

36.  35 

52 

59. 50 

38 

2.91 

1000 

9.75 

13 

26.  96 

15 

44.33 

39 

2.43 

1250 

7.44 

24 

21.71 

43 

35.  12 

88 

2.00 

1500 

6.32 

30 

17.  93 

10 

29.37 

42 

4.42 

25 


T  ABLE  V  i 


I 


s*1-. 

ending  Wave  Tub 

e  Measi 

:rsd  Da 

ta  -  Sssosat 

"t 

^  .  '  V  . 

£££ 3* 

First 

Min . 

Sec . 

Min. 

Third 

Min. 

Maximum 

HZ. 

cm. 

IV  . 

cm . 

mv. 

cm. 

m  v. 

mv. 

200 

4  1.6  1 

3.20 

51.5 

250 

31.74 

29.  5 

r ; 

s  ■ 

300 

26.20 

2.90 

83.  95 

3.  15 

55.5 

350 

21.51 

3.90 

70.04 

3.  30 

65.5 

400 

18.42 

10.  0 

61. 49 

9.40 

178 

450 

16. 1  7 

1.52 

54. 44 

1. 70 

26 

500 

14.26 

4.05 

48.  19 

3.40 

82.21 

4.  05 

61.5 

550 

1  1.10 

5.40 

42.  92 

5.30 

73.68 

5.60 

68.5 

600 

9.  18 

2.05 

37.  66 

1.98 

66.69 

1.70 

24.0 

650 

6.4  3 

6.05 

32.  79 

6.15 

60.39 

6.  10 

46.2 

700 

4.28 

4.00 

28.  62 

3.95 

52. 94 

3.80 

26.3 

750 

2.29 

2.90 

24.  54 

2.55 

47.49 

2.  64 

15.4 

800 

19.07 

3.95 

40.  41 

4.25 

62.37 

4.  15 

27.3 

850 

16.52 

1.10 

37.  32 

1.12 

57. 42 

1.  10 

10.3 

900 

14.41 

3.10 

33.  30 

2.90 

52.80 

3.  20 

47.5 

1000 

1  1.67 

.80 

29.  00 

.90 

46.18 

.  76 

22.5 

1250 

8.22 

.60 

21 . 86 

.45 

36.08 

.  52 

24.6 

1500 

6.53 

.65 

17.  96 

1.40 

29.49 

1.  13 

60.5 

After  the  wavelength  was  determined,  it  was  used 
to  scale  the  distance  to  the  first  node  (Table  VI)  .  These 
values  were  plotted  versus  the  frequency  (Figures  2. 4-2. 6). 
The  first  two  points  above  resonance  ware  also  plotted  as 

d/A  -  .5  (2.27) 

A  linear  regression  was  then  used  with  the  two  data  points 
both  above  and  below  resonance  to  determine  the  resonant 
frequency  "f  ".  The  resonant  frequency  for  each  resonator  is 
also  shown  on  Figures  2. 4-2. 6. 

Once  the  resonant  freguancy  was  obtained,  "XA" 
was  calculated  for  each  frequency  using  aquation  (2.13).  ”X, 

/ />  c"  was  then  calculated  using  equations  (2.24-2.26). 
Finally,  the  ratio  of  MXx//oc"  to  "X^”  (K^)  vas  calculated 
for  comparison  to  the  theoretical  value  (^)  determined  from 
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I  TABLE  71  ! 


Standing  Save  Tube  d/^  Data 


* 

I 

> 

ZL£2- 

Resonator  #1 

Resonator  #2 

Resonator  £3 

1 

1 

i 

Hz. 

d/* 

<3/A 

d/4 

I 

1 

i 

200 

.2416 

( 

250 

.2437 

.  2015 

.  2303 

1 

300 

.  220  1 

.  1726 

.2268 

1 

• 

350 

.21  34 

.  1028 

.2214 

I 

375 

.  0397 

1 

400 

.1934 

.  4683  (-.03  17) 

.2138 

I 

450 

.  1585 

.  3  86  5  (-  .  1 1  3  5) 

.2113 

1 

500 

.0954 

.  3454 

.2099 

| 

550 

.4848  (-.  0  152) 

.  1774 

I 

600 

.3906  (-.  1  094 

.  1596 

I 

650 

.3452 

.  1192 

1 

1 

700 

.3213 

.  0880 

1 

750 

.3094 

.  2943 

.0507 

1 

800 

.  4404(-.  0596) 

I 

850 

.  4039(-.  096  1 

I 

900 

.3754 

1000 

.2802 

.  2820 

.  3382 

1 

1250 

.2754 

.  2688 

.2950 

I 

| 

1500 

.2706 

.  2751 

.  2B44 

1 

9 

i _ _ _ _ _  _ _ » 


equation  (2.23).  Tables  VII-IX  list  the  experimental  values 
for  each  Helmholtz  resonator. 

c.  Error  Analysis 

Before  an  analysis  of  the  errors  could  be 
attempted,  it  was  necessary  to  identify  all  of  the  sources 
of  the  errors.  The  errors  were  a  result  of  the  measurements 
of  the  voltage  (Pm4A  ar.d  ,  frequency  and  distances  to  the 

nodes . 

Recalling  equation  (2.25) 

V,  =  cot h*'  (-&***---)  (2.25) 

» A 
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TABLE  71 I 


Standing  Ware  Tube  Reduced  Data  -  Resonator  “1  ! 


f,  =  542.7 

Hz  . 

Frequency 

c 

Zs£°£ 

X  1 
-m 

Hz. 

cm.  /s 

Pa-s/m 

s 1  X 

10“3 

250 

32050* 

-7.  39 

-928.1 

7. 96 

350 

34195 

-4.  15 

-491.5 

8.44 

400 

34432 

-2.63 

-336.3 

7.82 

450 

34605 

-1.  52 

-204.5 

7.43 

500 

34  80  0 

-.658 

-39.05 

7.  39 

550 

34430 

.094 

14.50 

6.48* 

600 

34  56  0 

.807 

109.  1 

7.40 

650 

34684 

1. 45 

196.9 

7.36 

700 

34636 

2.  07 

279.3 

7.41 

750 

3442  5 

2.55 

357.3 

7.  14 

1000 

34  80  0 

5.  21 

705.5 

7.  38 

1250 

34500 

6.22 

1014 

6.13 

1500 

34650 

7.  70 

1304 

5.90 

average 

34  56  0 

average 

7.  57 

deviation 

174 

deviation 

.39 

per  cent 

.50 

per  cent 

5.  15 

*  Values  not  used  for  the  calculation  of  the  average. 


it  can  be  seen  that  for  the  values  of  "P,^"  and  " " 
involved,  SHR  is  small  and  therefore,  an  error  in  the 

p 

standing  wave  ratio  (SWR  =  )  has  little  effect  on  the 

overall  calculations. 

Looking  at  equation  (2.24) 

_ _  _ _ _  (2.24) 

^  <i>,  +  TlS  Vl 

and  taking  into  account  that  "V,"  i3  small,  equation  (2.24) 
then  reduces  to 


_  AS-  &  -Jji'l  _ 
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TABLE  VIII 


Standing  Save  Tube  Reduced  Data  -  Resonator  * 2  \ 


387.8 

Hz . 

Freauer.cv 

c 

x  • 

—m 

Hz. 

cm. /s. 

Pa-s/m 

S~l  X 

10-3 

200 

34  39  2* 

-3.  86 

-551.9 

6.99 

250 

34  4  3  0  * 

-2.  43 

- 3o  1 .  o 

6.91 

300 

34  54  2 

-1.46 

-201.3 

7.25 

350 

344  4  7 

-.563 

-79.69 

7.06 

375 

33563* 

-.070 

-26.04 

2.69* 

400 

34  29  6 

.258 

24.03 

10. 74* 

450 

3351  6 

.745 

115.8 

6.43 

500 

34620 

1. 55 

199.2 

7.78 

550 

34557 

2.08 

276.6 

7.  52 

600 

34398 

2.  35 

349.4 

6.73 

650 

34437 

2.80 

418.6 

6.69 

700 

34  5  0  3 

3.  23 

485.2 

6.66 

750 

3450  8 

3.  72 

549.5 

6.77 

1000 

34  56  0 

5.  49 

84  9.6 

6.46 

1250 

34913 

8.51 

1130 

7.53* 

1500 

34470 

7.  11 

1400 

5.08* 

average 

3444  4 

average 

6. 94 

deviation 

314 

deviation 

.41 

per  cent 

.9  1 

per  cent 

5.9 

♦Values  not 

used  for 

the  calculat 

ion  of  the 

a  ver age 

or 


> 


-Ji i- 


/•* 


as  cot  Vi, 


away 

that 

from 

goes 

the 


from  resonance. 

On  the  basis  of  this  result,  it  was  determined 
the  best  data  was  obtained  in  the  regions  just  away 
resonance  for  the  following  reasons.  First,  since  tar. 
tc  zero  at  resonance,  "  terms  begin  to  dominate  and 
errors  in  "  VK 11  become  large  when  compared  to  " 


» 


> 


» 
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TABLE  IX  ! 

Standing  Wave  Tube  Reduced  Data  -  Resonator  #3  I 


f  „ =  773.0 

Hz. 

Frequency 

c 

i,o>z 

X  ’ 

—m 

Hz. 

cm.  /s. 

Pa-s/m 

S“ 1  X 

10-3 

200 

-7.86 

-2788 

2.  81 

250 

-6.93 

-2140 

3.  24 

300 

34650 

-5.  72 

-  1692 

3.  38 

350 

33999 

-5.  03 

-1357 

3.71 

400 

34  45  6 

-4.  09 

-1094 

3.74 

450 

34  44  3 

-3.  78 

-  878 

4.  31 

500 

33  97  5 

-3.63 

-  695 

5.22* 

550 

34  41  9 

-2.07 

-  536 

3.  86 

600 

34506 

-1.50 

-  396 

3.79 

650 

35074 

-.870 

-  269 

3.23 

700 

34062 

-.600 

-  154 

3.90 

750 

33900 

-.313 

-  46.7 

6.70* 

800 

3464  0 

.427 

53.  1 

8.04* 

850 

34  76  5 

.6  16 

147 

4.  19 

900 

34  55  1 

1. 06 

236 

4.49 

1000 

3451  0 

1.58 

402 

3.93 

1250 

34  82  5 

3.97 

772 

5.  14 

1500 

3444  0 

4.69 

1102 

4.26 

average 

344  5  1 

average 

3.74 

deviation 

327 

deviation 

.47 

per  cent 

.95 

per  cent 

12.6 

♦Values  not 

used  for 

the  calculation  of  the 

average 

Secondly,  as  the  frequency  gets  further  from  resonance  the 
cotangent  function  gets  larger  and  the  errors  also  get 
larger.  It  was  for  these  reasons  that  some  of  the  data  in 
Tables  VII-IX  were  deleted  for  the  calculation  of  the 
averages. 


This  completes 
wave  tube  technique.  The 
later  in  this  chapter  along 
dual-microphone  technique, 
technique  needs  addressing. 


the  discussion  of  the  standing 
final  results  will  be  examined 
with  those  for  the  Chung-Blaser 
One  additional  experimental 
The  resonant  frequency  of  each 


> 


» 
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Helmhcltz  resonator  was  measured  independently  of  e in  her 
method.  This  was  done  using  a  signal  generator  and  sc^hir 
to  excite  the  resonators.  Using  a  small  microphone  in  the 
resonators  ~o  detect  the  pressure  changes,  the  signal  gener¬ 
ator  was  tuned  to  the  frequency  which  resulted  in  the 
greatest  pressure  in  the  resonators.  These  resonant 
frequencies  ware  549.5  Hz.  for  resonator  1,  396.5  Hz.  for 

resonator  2  and  795.0  Hz.  for  resonator  3. 

C.  CHUNG-BLASER  DU AL-HICRO PHONE  TECHNIQUE 

The  dual- microphone  technigue  developed  by  Chung  and 
Blaser  provides  a  quick  method  for  measuring  the  acoustic 
impedance  and  reflection  coefficient  in  the  laboratory. 
However,  as  will  be  shown  later,  the  technique  suffers  seme 
degredation  in  accuracy  compared  with  the  results  obtained 
tv  the  SWT  method. 

1 .  Theory 

Chung  and  Elaser  have  shewn  [Ref.  1]  that  the 
complex  reflection  coefficients  at  the  two  microphone  loca¬ 
tions  can  be  expressed  in  terms  of  the  transfer  functions  of 
the  incident  and  reflected  waves.  This  expression  is 


J?.®  =  Atffl. 
HA) 


(2.28) 


where  (f)  and  R*  (f)  are  the  complex  reflection  coeffi¬ 
cients  evaluated  at  the  first  and  second  microphone  loca¬ 
tions,  respectively.  H*'(f)  and  Hr{f)  are  the  acoustic 
transfer  functions  associated  with  the  incident  and 
reflected  wave  components,  respectively,  evaluated  at  the 
two  microphone  loeatiens. 
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For  plane  wave  propagation  (no  mean  flew) 
neglecting  wall  losses,  H^(f)  and  Hr(c)  can  b*  exp;-''=-  a 

(2.  29) 


H<-  (f)  =  e'J*J 


and 


H r  (f)  =  e 


JkS 


(2. 30) 


where  "k"  is  the  wavenumber  and  "s'1  is  the  microphone 
spacing. 

Application  of  equation  (2.28)  results  in  the 
reflection  coefficient  "R"  on  the  surface  of  a  sample 
material  not  at  a  microphone  location  being  expressed  as 


=  .iLL.. 

R  V/ 


(2.31) 


whar  e 


H*  =  o'-**' 

l 


ar.d 


Hi  =  ejk' 


(2.32) 


(2.33) 


and  ”1"  is  the  distance  from  the  first  microphone  to  the 
surface  of  the  sample  material.  Osing  aquations  (2.31-2.33), 
the  reflection  coefficient  at  the  sample  is 


a  *  H,  eJzA/ 
where  [Ref.  1] 


(2.34) 


a  -  JtaS&rJi&L- 
•  '  //,<?)  - 


(2.35) 


H^f)  is  the  acoustic  transfer  function  between  microphones 
1  and  2  and  Ht(f)  and  H^(f)  are  the  transfer  functions  of 
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respectively,  between 


the  incident  and  reflected  waves, 
two  microphone  locations. 

Osirg  the  standard  relation  between  reflection 
i* c  —  e nt  end  cc  p .» e x  acoustic  tuipeda.*c^ 


<r  ^ 


-£i.  =  Ll+j£ L 
</-*) 


(2. 36) 


and  the  relationships  in  equations  (2.34)  and  (2.35),  the 
complex  specific  acoustic  impedance  can  be  determined  at  a 
distance  "l"  from  the  first  microphone  [Ref.  1] 


.it.  = .  I'jh&A  z5i  a 

/>C  3  CesCk  Km.  ^s(kl) 


(2.37) 


where 


-*£-  =  A.Ok) JvJj'tL-jflli -  (2  381 

j°c  Cos'CkC /-s3  <*  /-!)]}  *  /HJ^as'CV) 

H/t(f)  can  be  expressed  as 

H^(f)  =  (2.39) 


where  C/x  and  Q/J;_  are  the  real  and  imaginary  parts  of  the 
cross-power  spectrum  of  the  signals  at  microphones  1  and  2. 
S„  (f)  is  the  auto-power  spectrum  of  the  signals  at 
microphone  1. 

Finally,  Seybert  and  Ross  show  [Ref.  15]  that  the 
phase  change  upon  reflection  is  given  by 


0(f) 


(2.40) 


33 


where  Qtr  and  Ctr  ar€  the  imaginary  and  real  parts  of  t he 
cross-power  spectrum  of  the  incident  to  reflected  waves. 

2.  Equipment  Set-un 

Figure  2.7  shows  a  simplified  diagram  of  the 
Chung-Blaser  dual-microphone  technique  equipment  set-up.  At 
the  canter  of  the  system  are  the  Hewlett-Packard  HP-85 
computer  ar.d  HP-5420  Signal  Analyzer  System. 

a.  The  Hewlett-Packard  HP-85  Computer 

The  HP-85  computer  is  an  3  bit  microprocessor 
which  uses  the  BASIC  computer  language.  It  has  16K  bytes  of 
reai/write  memory  of  which  14,579  bytes  are  available  for 
use  by  the  operator.  The  unit  used  in  the  study  had  its 
memory  expanded  to  32K  bytes  by  a  memory  module.  The  unit 
has  a  mass  storage  capability  in  the  form  of  a  magnetic  tape 
reader/recorder.  It  also  includes  a  127  millimeter  diagonal 
black  and  white  electromagnetic  deflection  CRT  and  a  32 
character  per  line  thermal  printer/plotter.  To  interface 
with  other  equipment,  an  I/O  ROM  and  interface  card  were 
added  to  provide  HP-IE  (IEEE-488)  interface  capability. 

t.  The  Hewlett-Packard  HP-5423  Signal  Analyzer 

The  Hewlett-Packard  HP-5420A  Signal  Analyzer 
provides  digital  time-domain  and  frequency-domain  analysis 
of  lew  frequency  analog  signals.  The  HP-5420A  system 
includes  the  HP-5447B  Digital  Filter  Module,  HP-54410A 
Analog  no  Digital  Converter,  HP-5443A  Keyboard  and  Control 
Module  and  the  5441A  Display  Module.  The  equipment  used  in 
this  study  also  included  an  HP-10923A  HP-IB  interface  acces¬ 
sory.  This  allowed  its  use  with  an  HP-IB  controller.  It 
allowed  the  transfer  of  data  and  measurement  set-up  informa¬ 
tion  to  pass  to  and  from  the  controller  and  allowed  the 
reception  of  remote  control  commands.  Using  the  controller. 
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the  data  was  modified,  processed,  saved  and  recalled 
resulting  in  »  high  dears?  ci  flexibility. 

The  HP-5U20A  Signal  Analyzer  performes  a  variety 
of  rime-domain  and  frequency-domain  measurements.  ill  .‘■.in¬ 
urements  are  done  using  digital  signal  processing.  This 
involves  filtering  and  sampling  the  input  waveforms.  The 
analog  inpu*  is  first  filtered  by  a  lew-pass  filter  and  then 
passed  to  an  analog  tc  digital  converter.  Here,  the  waveform 
is  sampled  which  results  in  the  conversion  of  a  voltage  (at  a 
certain  point  in  time)  into  a  numerical  representation. 
These  numbers  are  then  filtered  and  processed  digitally  to 
produce  the  various  measurements  preformed  by  the  analyzer. 

c.  The  Impedance  Measurement  Tube 

The  tube  used  in  this  part  of  the  study  was  a 
piece  cf  3  inch  inside  diameter  PVC  pipe  12  inches  long  with 
a  3  inch  loudspeaker  mounted  at  one  end.  Two  holes  were 
drilled  in  the  tube  and  fitted  with  7/8  inch  Tygon  tubing  to 
hold  tha  two  microphones.  The  first  was  located  11.42  cm. 
from  the  sample  end  (opposite  the  loudspeaker)  and  the 
second  was  located  4.59  cm.  from  the  first.  The  distances 
were  measured  to  the  center  of  the  microphones. 

The  microphones  used  were  General  Radio  1962  1/2 
inch  elect ret-condenser  microphones.  These  microphones  have 
a  sensitivity  level  cf  -40  dB  re  1  v/Pa  and  have  typically 
flat  response  curves  from  20  Hz.  out  to  1800  Hz.  which  cover 
the  frequency  range  of  the  study. 

3 •  Experimental  Proced ures 

The  dual-microphone  technique  was  used  to  determine 
the  resonant  frequency  and  the  reactive  impedance  of  the 
same  three  Helmholtz  resonators  used  with  the  standing  wave 
tube  technique.  A  major  difference  between  the  two  techni¬ 
ques  made  itself  known  almost  immediately.  The  time  to 
collect  and  reduce  the  data  was  vastly  different. 
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Microphone  Calibration 


> 


a. 


Tae  two  aicrophone/preaaplifisr  channels  had  tc 
be  calibrated  for  phase  and  amplitude  differences  prior  to 
use  for  the  actual  measurements.  The  technique  used  for  this 
calibration  was  to  place  both  microphones  in  identical  sound 
fields  and  record  any  frequency-dependent  phase  and  ampli¬ 
tude  differences.  This  calibration  data  was  then  used  by  the 
computer  to  condition  the  data  prior  to  the  calculation  of 
reactive  impedance,  phase  angle  and  complex  reflection  coef¬ 
ficient.  It.  this  manner,  the  effect  of  differences  between 
the  two  microphone  channels  on  the  data  was  minimized. 
Additionally,  this  method  also  allowed  for  imbalances  in  the 
A/D  converter  channels  of  the  signal  processor. 

The  identical  sound  fields  were  obtained  by 
mounting  the  two  microphones  at  the  same  radial  position  at 
the  sample  end  of  the  tube.  This  was  based  on  the  assumption 
that  the  sound  field  is  radially  symmetric  at  frequencies 
balcw  the  sloshing  mode.  This  frequency  is  given  by 
CRef.  18] 


tu  *  «*c/d 

where  "  -s  -he  solution  to 

jCLd&ti  -  o 
d  ec 


(2.41) 


(2.42) 


and  "d”  is  the  diameter  of  the  tube.  The  sloshing  mode 
corresponds  to  which  is  equal  to  .5861.  This  results  in 
a  cut-off  frequency  of  2650  Hz.  The  calibration  program 
used  is  given  in  Appendix  A. 


b.  Validation  of  Tube  Performance 


the 

wall 


After  calibration  of  the  microphone  channels, 
microphones  were  placed  in  the  normal  positions  in  the 
of  the  tube  and  the  proper  operation  of  the  system  was 


» 


» 


i 


► 


» 


> 


» 


» 


» 


36 


verified,  A  4.03  inch  long  piece  of  ?VC  pipe  with  the  saa= 
dins  eter,  P  ro7i  <d  ed  vith  a  rigid  r.ila  z~.  i  ve  end,  wa-  '  - :  z  ~ 


a  sample  impedance.  The  expected  results  w  =  rs  a  rsfli 
coefficient  of  1  and  a  Known  phase  cna.nge  equal  no 


(2.«3) 


The  results  are  shown  in  Figures  2.3  and  2.5  Figure  2.8 
shows  that  the  reflection  coefficient  is  approximately  1  for 
frequencies  below  2.5  kHz.  except  at  1.5  kHz.,  were  the 
phase  goes  through  0  degrees  yielding  an  indeterminancv 
which  is  discussed  in  section  e.  which  follows.  Above  2.5 
kHz.  non-plane  wave  modes  degrade  the  results.  The  measured 
slope  of  the  phase  shift  versus  frequency  shown  in  Figure  X 
is  -53  degrees  +  1  degree.  Using  equation  (2.41),  this  slope 
produces  a  sound  speed  of  356  +  1.34  a/sec,  which  is  corsis- 
tant  but  significantly  higher  than  the  average  of  the  four 
values  obtained  by  the  S WT  technique  listed  in  Tables 
II,VII-IX  (c  —  345.  1  2  +  .38  m/sec). 

c.  Data  Collection 

Each  of  the  three  Helmholtz  resonators  was 
placed  at  the  sample  end  and  the  data  collected.  The  data 
was  tabulated  in  Tables  X-XII.  The  collection  and  recording 
of  the  data  was  done  completely  ay  the  computer  and  the 
program  is  given  in  Appendix  B. 

d.  Data  Reduction 

After  the  raw  data  was  collected  the  first  step 
was  tc  determine  the  resonant  frequency  of  the  Helmholtz 
resonator.  The  data  collected  allows  this  to  be  done  in  two 
different  ways.  Both  were  used.  First,  the  specific  reactive 
impedance  "Xs  // >c '•  was  plotted  versus  the  frequency  ar.d  a 
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TABLE  X 

Dual-Microphone  Measured  Data 


Resonator  #1  Rur.  1 


250 
350 
400 
450 
500 
550 
600 
6  50 
700 
750 
1  000 
1  250 
1  500 


Run  2 


Pa-s/m 


-5.  37 
-2.  83 
-1 . 91 
-  1.  23 
-.671* 
-.256* 
.  1  14* 
.395* 
.6  76 
.  9  19 
1  .86 
2.66 
3.03 


f.  =  587.  3  +-2.2  Hz. 


Phase  Angle 
deg. 


0 

23.9 
41.7 
65.  6 

105.  2 
147.  6 
166.  2 

135.6 

110.7 

93.9 
55.  2 
43.  2 
36.  5 


360  +  Phase  kr\o. 

deg. 


360 
336.  1 

318. 3 

294. 4 
254. 8 

212.4 


f#  =  590.6  t  *2  HZ. 


ls6>z 

Ph^se  Angle 

36  0  ♦  Phase  Ancj 

HZ. 

Pa-s/m 

deg. 

deg. 

250 

2.  13 

19.  2 

379.  2 

300 

-1 . 58 

-  8.  4 

35  1 . 6 

350 

-2.  14 

-  27.4 

332 . 6 

400 

-  1 . 64 

-  52.2 

307.  8 

450 

-1. 07 

-  78.6 

281. 4 

500 

-.591 

-113. 5 

24  6 .  5 

550 

-.2  24 

-152.  6 

207. 4 

600 

.076 

170.  8 

650 

.3  50 

140.  3 

700 

.6  29 

1 14.  6 

750 

.8  90 

9  6.0 

8  00 

1 .  19 

79.9 

900 

1.  52 

6  6.5 

1000 

2.22 

4  7.  8 

1250 

2.  96 

37.  1 

1  500 

3.  17 

34.9 

linear  regression  was  used  done  using  the  two  points  or. 
either  side  of  the  origin  (4  points  total)  where  the  r--~- 
tiva  impedance  was  equal  to  zero.  This  gives  the  resonant 
frequency  sinca  the  reactive  impedance  of  a  resonator  equals 
zero  at  resonance.  The  graphs  ana  resonant  frequencies  for 
each  resonator  are  shewn  in  Figures  2.10-2.14 

The  second  method  for  determining  the  resonant 
frequency  involved  the  use  of  the  phase  angle  between  the 
incident  and  reflected  waves.  When  the  phase  angle  was 
plotted  versus  frequency  (Figures  2.15-2.19)  resonance  was 
obtained  from  the  graph  at  the  point  where  the  phase  angle 
was  180  degrees.  Again  a  linear  regression  was  used  with  the 
2  points  on  either  side  cf  the  180  degree  phase  angle. 
Figures  2.15-2.19  also  show  the  frequency  determined.  The 
resonant  frequency  used  in  subsequent  calculations  was  the 
average  of  the  two  frequencies  determined  from  the  two 
methods.  The  reproducibility  cf  the  results  can  be  estimated 
from  the  ratio  of  the  standard  deviation  to  the  average  for 
the  four  (two)  separate  measurements  of  the  resonant 
frequencies  for  resonators  1  and  2,  which  were  0.5?  and 
2.0%,  respectively. 

Once  the  resonant  frequency  was  determined,  it 
was  used  along  with  equation  (2.13) 


—  £  _ 

—  i. 


C/f 


(2. 13) 


to  obtain  "X^".  Using  this  along  with  the  reactive  impedance 
ratio  "Xs//>c»,  the  experimental  value  for  the  ratio 


*w> 


>C 


(2.42) 
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!  TABLE  XI 

i  Duai-aicrophone  Measured  Data  -  aesonacor  #2  Run  1 


fa  =  435.8  +  3.8 

Hz. 

lL23- 

phase  Ar.cjie 

360 

+  Phase  A  r.  c 

Hz. 

Pa-s/m 

deg. 

deg. 

20  0 

.215 

7.5 

367. 5 

250 

-4. 24 

-  1.6 

358. 4 

300 

-.151 

-  8.  3 

351. 7 

350 

-.4  36 

-  47.3 

312. 7 

4  00 

-.2  06 

-121.7 

238. 3 

450 

.  1  00 

164.  6 

500 

.3  86 

133.  0 

550 

.667 

109.7 

600 

.  8  93 

94.  8 

650 

1.07 

85.0 

700 

1.  27 

76.0 

750 

1  .  40 

70.4 

800 

1.  60 

63.8 

900 

1  .  84 

56.7 

1000 

2.  12 

50.  1 

1  250 

2. 57 

42.2 

1  500 

2. 83 

38.  1 

Run  2 

f . *  421.7  * 

.  3  Hz. 

£isa* 

Phase  Angle 

360 

+  Phase  Ang 

HZ. 

Pa-s/m 

deg. 

• 

tj* 

0) 

200 

.  163 

1.6 

361 . 6 

250 

-2. 03 

-  41.7 

318. 3 

300 

-.987 

-  87.4 

272.6 

350 

-.4  34 

-132.  5 

227.  5 

400 

-.0  82 

-170.5 

189. 5 

450 

.171 

160.0 

500 

.  4  12 

134.4 

550 

.5  93 

1  16.  9 

600 

.  8  13 

100.7 

650 

.993 

89.8 

700 

1.  27 

75.7 

750 

1. 50 

67.  1 

800 

1 .  78 

58.  4 

900 

1  .  91 

54.8 

10  00 

2.65 

41.2 

1  250 

3.35 

33.  2 

1  500 

3.38 

33.0 

40 


ft 


TABLE  XII 

Ou&l-ii  icrophcne  Measured  Data  -  Sescnatnr  *3 


f,  =  947. 1  *  1.0  Hz. 


Freq. 

Phase  Ancle 

360  +  Phase 

Hz. 

Pa -s/s 

deg. 

deg. 

400 

.  7  19 

22.  1 

38  2.  1 

450 

.  1  74 

5.6 

365 . 6 

500 

-.0  39 

.  9 

359 .  1 

550 

-.5  97 

-  1  5.  4 

344 . 6 

600 

-1.01 

-  25.4 

334 . 6 

650 

-1  .  06 

-  34.6 

325. 4 

700 

-.921 

-  56.4 

303 .  6 

750 

-.8  34 

-  69.4 

290 . 6 

800 

-.6  25 

-  94.6 

265. 4 

850 

-.3  80 

-126.9 

233 .  i 

900 

-.201 

-150.7 

209 .  3 

950 

.048 

173.  5 

1000 

.  1  96 

155.  2 

1  050 

.3  85 

135.5 

1  100 

.571 

113.  1 

1200 

1 .05 

8  6.5 

1300 

1.34 

7  3.0 

1400 

1 . 57 

64.7 

1500 

1 .  69 

6  1.2 

ft 


l 


» 


I 


ft 


was  calculated  for  comparison  with  the  theoretical  value 

obtained  from  equation  (2.23)  ft 

K  _  (2.23) 

C/ff 

ft 

Tables  XIII-XV  show  the  results  foe  "K^p"  for  each  of  the 
resonators  tested. 

e.  Error  Analysis 

The  Chung-Biaser  dual-microphone  technique  had 
some  inherent  errors  built  into  it  by  virtue  of  the  equip¬ 
ment  construction  and  mathematical  techniques.  First,  the 

» 

4  1 


ft 


TABLE  XIII 


I  Dual-Microphone  Beducad  Data  -  Resonator  #1  Sun  1 


Frequency 

HZ. 

s~  i 

x  10- 

250 

-1142 

5.14* 

350 

-  644.2 

4.  39 

400 

-  470.0 

4.06 

450 

-  323.3 

3.80 

500 

-  196.0 

3.42 

550 

-  82.7 

3.  10* 

600 

20.0 

5.69* 

650 

1  14. 6 

3.45 

700 

202.9 

3.33 

750 

286.0 

3.21 

1000 

652.0 

2.85 

1250 

971.6 

2.74 

1500 

1268 

2.39 

average 

3. 52 

deviation 

.79 

per  cent 

22.4 

Run  2 

Frequency 

li 

HZ. 

S“  1 

x  10- 

250 

-1146 

1.86* 

300 

-  863.1 

1.83* 

350 

-  646.9 

3.31 

400 

-  472.3 

3.  47 

450 

-  325.4 

3.29 

500 

-  197.9 

2.97 

550 

-  84.4 

2.65* 

600 

18.5 

4.  12* 

650 

113.2 

3.09 

700 

201. 5 

3.  12 

750 

284. 8 

3.  13 

800 

363.8 

3.  27 

900 

512.3 

2.97 

1000 

651.  1 

3.41 

1200 

909.2 

3.  26 

1400 

1151 

2.75 

average 

3.  17 

deviation 

.21 

par  cent 

6.62 

*Values  not  used  for  the  calculation  cf  the  averag 


TABLE  XIV 

Dual-Micro? bone  Seduced  Data 


Resonator  *2  Pur.  1 


requency 

Hz  . 

200 

250 

?  *>  f% 

35  6 
400 
450 
500 
550 
6  00 
650 
700 
750 
30  0 
900 
1000 
1200 
1400 


Run  2 


200 

250 

300 

•350 

400 

450 

500 

550 

600 

650 

700 

750 

800 

900 

1000 

1200 

1400 


u 

iW 

S“  1 

x  1 0-  3 

-  749.6 

.  29* 

-  509.7 

8.  32* 

-  33  3.1 

.  45* 

-  192.6 

2.  26 

-  74. 8 

2.75 

28.  0 

3.58* 

120.2 

3.21 

204.  7 

3 .  26 

283.  5 

3.  15 

357.  8 

2.99 

428.7 

2.  96 

496.8 

2. 82 

562.  6 

2.84 

689.0 

2.67 

810.1 

2.62 

1042 

2.47 

1264 

2.  24 

average 

2.79 

deviation 

.  33 

per  cent 

11.8 

X* 

—  «-»/0 

s~ 1 

x  10-3 

-  689.6 

.  24* 

-  46  1.7 

4.40* 

-  293.1 

3.37 

-  158.3 

2.74 

-  44.  8 

1 . 83* 

54.6 

3.  13 

144.  2 

2. 86 

226.  5 

2.62 

303.  5 

2.  69 

376.  3 

2.64 

445.  8 

2.85 

512.8 

2.93 

577.6 

3.08 

702.3 

2.72 

822.  1 

3.  22 

1052 

3.  18 

1273 

2.66 

av erage 

2.91 

deviation 

.25 

par  cent 

8.59 

used  in  the  calculation  of  the  averag' 


♦Values  no 


r 


TABLE  XV  , 

! 

Du ai-Hicro phene  Reduced  Data  -  Resonator  *3  i 

! 


Frequency 

'^e*P 

Hz. 

S“  1 

x  10-3 

400 

-  1844 

.  39 * 

450 

-  1545 

.11* 

500 

-1295 

.03* 

550 

-  1082 

.  06* 

600 

-  895.9 

1.  13 

650 

-  730.9 

1.45 

700 

-  582.2 

1 . 58 

750 

-  446.8 

1.87 

800 

-  322.0 

1.94 

850 

-  206.0 

1.84 

900 

-  97.3 

2.07 

950 

5.  2 

9.24* 

1000 

102.4 

1.91 

1050 

195.  2 

1.97 

1100 

284.0 

2.01 

1200 

452.0 

2.32 

1300 

609.  6 

2.  20 

1400 

758.  9 

2.07 

1500 

901.6 

1 . 87 

average 

1.91 

deviation 

.  23 

per  cent 

12.0 

*  Values  not  used  in  the  calculation  of  the  average. 
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construction  of  the  tube  limits  the  range  of  frequencies 
over  which  good  data  could  be  obtained.  The  diameter  of  the 
tube  determined  the  upper  limit  of  measurable  frequency. 
Since  the  tube  was  valid  only  for  a  plane  wave,  the  frequen¬ 
cies  measured  had  tc  be  below  the  curoff  frequency  for  the 
tube.  For  the  tube  used  in  this  study  the  cutoff  frequency 
was  approximately  2600  Hz.  The  construction  of  the  tube  also 
limits  the  effectiveness  of  measurements  at  the  lower 
frequencies.  The  critical  parameter  in  this  case  is  the 
spacing  between  the  microphones.  For  a  given  separation  as 


44 


the  wavelength  aets  longer  (frequency  lower),  the  differ¬ 
ences  nr  the  phase  and  amplitudes  at  the  two  inr  ?  r'~  ~n- 
locatior.s  became  small  ar.d  the  noise  dominated.  The  lov 
frequency  limit  for  the  UDe  was  about  200  Hz. 

Finally,  the  mathematical  technique  resulted  in 
an  indeterminacy  in  the  reflection  coefficient.  This 
occurred  when  the  phase  angle  was  equal  to  zero.  when  this 
occurred  equation  (2.35)  reduced  to 

S  =  0/0 

wnich  is  indeter minant.  The  result  was  that  data  obtained 
for  the  frequencies  about  50  Hz  on  either  side  of  the  point 
where  the  phase  angle  equalled  zero  was  greatly  in  error. 

D.  COMPARISON  OF  THE  TWO  TECHNIQUES 

The  two  techniques  were  consistant  with  each  ether  and 
with  the  expected  theoretical  values.  First,  the  techniques 
will  be  locked  at  for  comparison  with  the  the  theoretical 
results  and  then  for  comparison  with  each  other. 

1 .  Standing  Wave  Tube  Technique 

Figures  2.20-2.22  show  the  comparison  of  the  experi¬ 
mental'  and  theoretical  values  for  "X,//>c".  The  theoretical 
values  were  obtained  using  "Ky*"  (eqn.  2.23)  and  the  calcu¬ 
lated  values  for  "X,*  ”  (Tables  VII-IX)  .  These  were  tabulated 
in  Tables  VII-IX  As  can  be  seen  in  Figures  2.20-2.22  the 
values  compare  favorably  with  a  divergence  noted  at  the 
higher  frequencies  away  from  resonance.  This  divergence  was 
explained  in  the  error  analysis  in  section  B.3.c. 

Table  XVI  shews  the  resonant  frequency  and  both 
constants  for  each  resonator,  along  with  the  resonant 
frequency  determined  by  the  independent  method  discussed  in 
section  B.3.c.  As  can  be  seen,  and  were  within 
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the  calculated  error  for  each.  The  average  per  cent  error 
between  the  tvc  values  for  "K  "  was  found  to  be  4.1?.'  for 
the  standing  wave  tube  technique.  The  average  per  cer.-  -rrcr 
for  the  resonant  frequency  was  found  to  be  -2.23%. 

2.  Chung- Blase r  Dual-M icrophone  Technique 

Figures  2.23-2.27  show  the  comparison  of  the  experi¬ 
mental  and  theoretical  values  for  11 Xs//?c ".  The  theoretical 
values  were  obtained  in  the  same  manner  as  in  the  SWT  tech¬ 
nique.  The  qraphs  shew  a  good  agreement  except  at  the  lower 
and  higher  frequencies  away  frem  resonance  for  reasons 
previously  noted  in  section  C.3.e. 

Table  XVI  shews  the  resonance  frequency  and  beth 
values  obtained  for  the  dual-microphone  technique.  The 
values  determined  fer  and  "K^  ••  do  not  in  all  cases 

overlap  within  the  areas  of  uncertainty  as  was  the  case  with 
the  SWT  technique.  As  a  result,  the  average  per  cent  error 
((Kw<(  -  K^)/K^)  x  100  was  found  to  be  18.4%.  The  average  per 
cent  error  between  the  independently  measured  resonant 
frequency  and  that  measured  by  the  dual-microphone  method 
was  found  to  be  +11.4%. 

3.  Standing  Wave  Tube  vs.  Du^l- Microphone 

The  greatest  discrepancy  between  the  two  methods  was 
the  measured  resonant  frequencies.  The  resonant  frequency 
measured  by  the  dual-microphone  technique  was  always  higher 
than  the  resonant  frequencies  measured  by  the  other  two 
methods.  The  average  error  of  +11.4%  was  significantly 
larger  than  the  error  for  the  standing  wave  tube  method  of 
-2.23%.  Additionally,  the  determination  of  the  constants 
for  the  dual-microphcne  method  produced  larger  discrepancies 
between  the  theoretical  and  experimental  values.  Finally, 
the  calculated  sound  speed  for  the  dual-microphone  method 
was  3.  235  higher  than  the  SWT  method  which  was  in  good 
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TABLE  XVI 

Ccmoarison  of  Results 


HZ. 

sac. 

x  10-3  X 

sac. 

10-3 

% 

z. 

54  2. 7 

7.35  +  .18 

7.57 

♦  .39 

♦  2.99 

589.  5 

3.67  t  .10 

3.52 

±  *79 

♦  4.09 

590. 5 

3.66  +  .10 

3.17 

i  *21 

-13.4 

Freq , 
HZ. 

RES  jM 

549.5  £  1.5  HZ. 

S  XT  54  2 

Dual  run  1  58  9 

Dual  run  2  590 

RES  #2 

398.5  +  3.5  HZ. 

SWT  38  7 

Dual  run  1  432 

Dual  run  2  422 

RES  #3 

795.0  +  4.0  Hz. 
SWT  77  3 


PER  CENT  ERROR  IN  RESONANT  FREQUENCY 

I  Re  gen  at  or  JM  -  549.5  Hz. 

SWT  -  542. 7  HZ .  -  1.24* 

Dual  -  590.3  Hz.  +7.42% 

Resonator  #2  -  398.5  Hz. 

SWT  -  387.4  Hz.  -2.79* 

|  Dual  -  428.8  Hz.  +7.60* 

!  Resonator  #3  z  221*5  25* 

SWT  -  773. 0  Hz .  -  2.77% 

Dual  -  974. 4  Hz .  ♦  22.6% 


x  1001 


387.  4 

7.37 

±  *15 

6.94 

♦ 

.41 

-5.83  | 

432.  0 

3.4  5 

±  .33 

2.79 

+ 

.33 

-19.1 

422.  0 

3.6  8 

.07 

2.91 

1 

.  25 

-20.9  1 

I 

1 

1 

Iz.  1 

1 

77  3.  0 

3.6  1 

t  .09 

3.74 

.47 

♦  3.60  1 

) 

974. 4 

1.4  2 

±  .04 

1.91 

* 

.23 

♦34.5  I 

1 

Figure  2.3  Standing  Have  Tube  Equipment  Set-Op. 


Figure  2.4  Resonator  #1  Graph  of  &/$  vs.  Frequency 


Figure  2.5  Besonator  #2  Graph  of  d/q  vs.  Freguenc 


PRE-AMP  - »  HP-54  4 10/1 - J  HP-5420  Hp_85  lIp_ 


Figure  2.7  Dual-Microphone  Technique  Equipment  Set-Up. 
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HLS //>c"  vs.  Freq.  -  Resonator  #1  Run 


Figure  2.11  Dual-tlicrophone  "Xs  //oc"  vs.  Freg.  -  Resonator  #1  run 


icrophone  uXs//oc"  vs.  Freg.  -  Resonator  #2  Run 


Pigure 


313NV  3SVHd 


Figure  2.16  Dual-Microphone  Phase  Ang.  vs.  Freg.  -  Resonator  #1  Run 


•  • 


65 


Figure  2.17  Dual-Microphone  Phase  Ang.  vs.  Freg.  -  Resonator  #2  Run 


67 


Figure  2.19  Dual-Microphone  Phase  Ang.  vs.  Freq.  -  Resonator  #3. 


Comparison  of  Reactive  Impedances  -  Resonator  #1. 


Figure  2.21  SWT  Comparison  of  Reactive  Impedances  -  Resonator  #2. 


Figure  2.23  Dual-Mic.  Comparison  of  Reactive  Impedances 


Figure  2.2<»  Dual-flic,  comparison  of  Reactive  Impedances 


Figure  2.25  Dual-Mic.  Comparison  of  Reactive  Impedances 


H- 

Z 


Pigure  2.27  Dual-Hic.  Comparison  af  Reactive  Impedances 


III.  CONCLUSIONS 


Throughout  this  study  both  techniques  displayed  advan¬ 
tages  and  disadvantages.  The  criteria  for  saying  one  is 
"better"  than  the  other  depends  upon  the  purpose  fcr  which  # 

the  information  is  needed.  No  attempt  will  be  made  to  decide 
which  technique  is  "better".  The  strengths  and  weaknesses 
will  be  pointed  out  and  the  potential  users  of  this  informa¬ 
tion  can  decide  which  technique  more  closely  fits  their  , 

needs. 

It  was  clear  from  the  results  in  Table  XVI  that  the 
standing  wave  tube  technique  was  the  more  accurate  when 
compared  tc  the  expected  theoretical  values  for  the  » 

Helmholtz  "standards".  In  addition,  the  SWT  technique  had  a 
higher  average  relative  precision  than  the  dual-micrcphcne 
technique.  The  average  fcr  the  SWT  was  4.2%  whil“  the 

average  relative  precision  for  the  dual-microphone  technique  § 

was  18.4%.  We  believe  the  reason  for  the  factor  of 
fcur-and-one-half  better  results  for  the  SWT  technique  was 
due  to  one  major  characteristic  of  the  two  techniques.  The 
standing  wave  established  in  the  two  tubes  were  thsoreti-  * 

cally  the  same.  However,  the  locations  where  the  appropriate 
measurements  were  obtained  were  different  for  each  tech¬ 
nique.  Fcr  the  SWT  technique  the  measurements  were  always 
taken  at  the  optimum  locations,  ie.  a  maximum  or  a  minimum  ■ 

pressure  point.  Therefore,  it  was  possible  to  obtain  the 
best  data  available.  For  the  dual-microphone  technique  the 
measurements  were  taken  at  locations  on  the  standing  wave 
predetermined  by  the  microphone  locations.  These  locations  • 

were  only  optimized  when  the  microphone  separation  was  equal 
to  an  odd  multiple  of  a  qu arter- wavelangt h  of  the  frequen¬ 
cies  measured.  Therefore,  the  data  obtained  was  not 
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necessarily  the  best  avails  ole  and  the  result  was  a  ra¬ 
tion  in  accuracy  and  relative  precision  for  the  :aal- 
nicrcphbne  technique.  This  does  not  provide  an  obvious 
explanation  for  the  error  in  measured  sound  speed  or  signif¬ 
icant  error  in  measured  Helmholtz  frequency  obtained  by  the 
dual-microphone  method. 

It  would  appear  that  the  SWT  method  was  the  more  desir¬ 
able  based  on  accuracy  and  relative  precision.  This  may  nor 
be  the  case.  When  numerous  samples  were  involved  or  the 
results  desired  over  a  large  frequency  range,  the  SWT  method 
became  tedious  and  cumbersome.  It  involved  the  manual  taking 
and  reduction  of  data.  The  dual-microphone  technique  was  a 
great  deal  faster  and  easier  to  use.  The  time  to  obtain  the 
results  for  one  sample  over  a  frequency  range  of  200-1500 
Hz.  was  approximately  20  minutes  for  the  dual-microphcne 
technique  and  6  hours  for  the  SWT  technique. 

Thus  when  choosing  which  technique  is  more  desirable  for 
a  given  application,  it  must  be  decided  whether  speed  or 
accuracy  is  more  important.  The  degradation  in  the  accuracy 
and  relative  precision  for  the  dual-microphone  technique 
might  be  considered  minor  in  those  applications  where  speed 
is  important. 
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IV.  HECOKSEN  DAT  loss  FOH  I«Pf07EJ!5NTS  AND  FU5THEP 

While  the  experiments  conducted  in  this  study  were 
straightforward,  several  changes  could  be  made  to  improve 
both  techniques.  The  biggest  change  would  be  to  adapt  the 
SWT  technique  for  computer  control.  This  could  be  accom¬ 
plished  using  a  servo-mechanism  with  a  worm-drive  on  the 
microphone  car.  The  computer  could  then  automatically  set 
the  frequency,  determine  the  maximum  and  minimum  values  and 
their  locations  and  compute  the  necessary  information.  This 
would  greatly  enhance  the  SWT  technique. 

One  improvement  or  change  for  the  dual-microphone  tech¬ 
nique  that  needs  investigation  is  the  use  of  a  pure  tone  as 
the  signal  instead  of  random  noise.  It  would  be  possible  to 
step  through  the  frequencies  desired.  The  use  of  a  pure  tone 
would  allow  the  use  of  a  narrow  band  filter  and  could  result 
in  a  reduction  in  the  noise  interference  and  would  include  a 
useful  redundancy  as  the  wavelength  and  thus  phase  differed 
between  the  two  microphones. 

Although  hare  we  only  dealt  with  reactive  impedance,  it 
would  also  be  possible  to  extract  an  absolute  value  of  the 
resistive  component  of  the  Helmholtz  resonator  with  only  one 
additional  measurement.  If  the  quality  factor  of  the  reso¬ 
nator  is  measured  and  the  absolute  reactive  impedance  is 
known,  the  absolute  resistive  impedance  is  fixed.  This  value 
can  then  be  compared  with  the  impedance  measurement. 

Finally,  several  changes  could  be  made  to  adapt  the 
dual-microphone  technique  for  use  in  a  water  filled  tube. 
First,  the  tube  should  be  constructed  with  numerous  hydro¬ 
phone  locations  so  that  the  separation  between  hydrophones 
could  be  varied  somewhat  to  optimize  the  data  for  a  given 
frequency  range.  Additionally,  the  use  of  a  bubble  of  known 
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dimensions  could  be  used  ;o  determine  the  accuracy 
tive  precision  of  the  water  tube  system.  This  would 
ogous  to  the  Helmholtz  resonator  used  in  this  study 
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4£PEN°IX  A 

MICROPHONE  CALIBRATION  PROGRAM. 

1000  cpticn  base  1 

1010  !  matrix  'a*  holds  HP5420  data  block  header 

1020  !  matrix  'A'  holds  auto  ar.d  cross  spectrum  data 
1030  !  matrix  '  C'  is  calibration  vector 
1040  !  address  of  HP-5420  is  704 

1050  dim  W  (2,  16)  ,  A  (3,256)  ,  C  (2 ,256) 

1060  output  704  ;  "R  S"  3  S=SP0LL  (704)  3  wait.  500  3  S#0  the 
1060 

1070  !  HP5420  test  setup: 

1080  !  60  stable  averages 

1090  !  channel  1  active 

1100  !  transfer  function 

1110  !  randon  signal 

1120  !  1.6kHz  bandwidth 

1130  output  704  ; " 6 0 , 1  A V  CH  5  FR  2  SG  1.6  KV  3W  0  CH  0  AC 
10  RG" 

1140  S  =  SPOLL  (704)  3  wait  500  3  if  S  *)  then  1140 

1150  if  S#0  then  1 1  40 

1160  !  start  data  collection 

1170  disp  "press  end  line  to  continue"  3  input  AS 
1180  output  704  ; " ST" 

1190  S=SPOLL  (7 04)  3  wait  500  3  if  S#63  then  1190 
1290  disp  "press  end  line  to  continue"  3  input  AS 
1210  gosub  1330  3  gosub  1600 

1220  !  create  data  file'CALDAT'  if  not  present  on  current 

tape 

1230  on  error  goto  1240  3  create  "CALDAT" , 1 , 4096 
1240  off  error  3  assign#  1  to  "CALDAT" 

1259  print  *  1,1  ;  C  (,) 

1260  I  close  file  to  record  data  on  tape 


127  0  assign#  1  to  * 

1280  disp  "calibration  data  snored" 

1290  disp  "end  program" 

1300  clear  7 
1310  stop 
1320  and 

1330  !  collect  data 

1340  !  transfer  aatc  spectrum  data  ::  KP85 

1350  output  704  ; " 1  TC  1  PM  1  A0" 

1360  S=SPOLL  (704)  2  wait  500  2  if  S#4  then  1360 

1370  output  704  ; "5  01SA" 

1380  S=SPOLL  (704)  2  wait  500  2  if  S#32  then  1380 

1390  for  1=1  to  16 

1400  !  enter  5420  header  information  into  matrix  w 

14  10  enter  704  ;  H(1,l) 

1420  next  I 

1430  !  enter  data  in  matrix  A 

1440  for  1  =  1  to  256  a  enter  704  ;  A  (1,1)  2  next  I 
1450  !  transfer  cross  spectrum  data  to  HP85 

1460  output  704  ;»RS" 

1470  S=SFOLL  (704)  2  wait  500  2  if  S#32  then  1470 

1480  output  704  ;"CP." 

1490  S  =  SPOLL  (704)  2  wait  500  2  if  S#0  then  1490 

1500  output  704  ;"501SA" 

1510  S  =  SPOLL  (704)  2  wait  500  2  if  5*32  then  1510 

1520  for  1=1  to  16  2  enter  704  ;  tf(2,I)  2  next  I 

1530  for  1=1  to  256 

1540  !  enter  real  cross  spectrum  data  to  matrix  'A' 

1550  enter  704  ;  A (2,1) 

1560  !  enter  imaginary  cross  spectrum  data  to  matrix  'A' 

1570  enter  704  ;  A  (3,1) 

1580  next  I 
1590  return 

1600  !  calculate  calibration  vector 

1610  print  "calculating" 
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DO  AL “MICROPRO  HE  COMPUTER  PROGRAM 


10 

laddress  of 

HP5420  si  704,  plotter  is  705 

20 

!  list  of  variables 

30 

'■  A  ( 1 ,) 

in  ic  1  a  a 

40 

!  A  (2  , ) 

cross  spectrum,  real  part 

50 

!  A  ( 3  , ) 

cress  spectrum,  imaginary  part 

6  0 

!  C  ( 1  r) 

calibration  vector,  real  part 

70 

!  C<2,) 

calibration  vector,  imaginary  part 

80 

!  RH 

real  part  of  acoustical  transfer  function 

90 

!  IH 

imaginary  part  of  acoustical  transfer 

fa net  ion 

100 

!  A1() 

reflection  coefficient 

110 

!  PHI  () 

phase  change 

120 

!  B  O 

impedance,  real  part 

130 

l  &() 

impedance,  imaginary  part 

140 

option  base 

1  3)  dim  A  (3,256)  ,A  1  (256)  ,C(2,256)  ,M(256), 

R  (256)  ,PHI  (256) 

150 

l  all  following  air  constants  and  calculations  are  in  SI 

units 

160 

T  3=29  5.15 

!  temperature 

170 

G  3=  1 . 4 

l  specific  heat  ratio 

180 

P3  =  1 0  132  5 

!  air  pressure 

190 

R3=287.0537 

!  gas  constant 

200 

C=SQR  (G3*R3*T3)  !  speed  of  sound  in  air 

210 

RO-P3/  (R3*T3)  !  air  density 

220 

ROC=RO*C 

230 

assign#  1  to 

"CALDAT"  d)  read#  1,1  ;  C(,) 

240 

disp  "enter 

number  of  materials  to  be  measured" ; din  put 

13 

250 

for  19=1  to 

18  a)  print  "we  will  now  measure  material  #" 

:  19 

•  • 


»  •- 


»  • 
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260  disp  "if  sample  is  ready,  press  end  line"  3  ir.cn-;  A  J 

270  cut  par  704  ;  "RS"  a  S=SPQLL  (704)  3  wait  500  3  if  5  SO 

-her.  270 

230  output  704  ;"30,1  A  7  1  CK  5  FS  2  SG  1.6  KV  B.i  0  CH  0  AC 
AC  10  RG" 

290  S=SPCLL  (704)  2)  wait  500  5  if  S*0  then  290 

300  output  704  ;»ST  MR" 

310  S=S  FOIL  (704)  3  wait  500  3  if  S#68  then  310 
320  disp  "and  data  collection" 

330  output  704  ; ” 1  FM  1  TC  1  AU" 

340  S=S  POLL  (704)  3  wait  500  3  if  S#4  then  340 

350  J=1  3  gosub  530 

360  output  704  ; "RS"  3  S=SPOLL  (704)  3  wait  500  3  if  S#0 
then  360 

370  output  704  ; "CR" 

380  S=SPOLL  (704)  3  wait  500  3  if  S#0  then  380 

390  gosub  530 

400  output  7 04  * ;  "RS"  3  S=SPOLL  (704)  3  wait  500  3  if  S#0 
then  400 

410  disp  "off  to  calculate!  this  will  take  several 
minutes" 

420  gosub  600 

430  disp  "press  end  line  to  store  impedance  data  on  tape" 
440  input  A  $ 

450  !  the  following  statements  save  data  on  tape 

460  create  " IM  PD  AT  ",  256, 2  4 

470  assign#  1  to  "IMPDAT"  3  for  1=1  to  256  3  F=(I-1)*D1  3 
print  #1, 1 ;  F,  R  (I)  3  next  I  3  assign#  1  to  "IMPDAT2" 

480  disp  "end  data  collection  for  set" ; 19 

490  next  19 

500  disp  "end  program" 

510  end 

520  !  start  data  collection 

530  output  704  ; "5  0  ISA" 

540  S=SPOLL  (704)  3  wait  500  3  if  S#32  then  540 
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550  for  1=1  to  16  2  enter  704  ;  w  a  next  I 
560  if  J=2  then  580 

57  0  for  1  =  1  to  256  3  enter  704  ;  A (J,I)  a  next  T  2  crlrt 
"auto-spectral  data  t ransf erred"  3  J=J+1  3  return 
580  for  1=1  to  256  3  enter  704  ;  A  (2,1)  3  enter  704  ;  A  (3,1) 
3next  I  3  print  "cross-spectral  data  transf  erred"3return 
590  !  calculate  impedance 

60  0  disp  "calculating" 

610  for  1=1  to  256 

620  C1=C(1  ,I)*C  (1,1)  +C  (2,1)  *C  (2,1) 

630  TA=  (A  (2,  I)  *C  (1,1) + A  (3,1)  *C(2,I)  )  /Cl 
64  0  A  (3,1)  =  ( A (3,  I)  *C  (1,1) -A  (2, 1)  *C(2,I))/C1 
650  A  (2 , 1)  =T  A 

660  next  I 

670  !  L  is  the  distance  from  the  sample  to  the  center  of 

the  farthest  mic. 

680  !  S  is  the  spacing  between  the  microphone  centers 

690  !  D1  is  the  frequency  spacing  of  the  5420  samples 

700  S=.0459 

710  L= . 1 1 42  3  D1 =6  .25 

720  !  calculate  transfer  function,  impedance,  reflection 

coefficient  and  phase  change 
730  for  1=1  to  256 
740  F=  (I- 1 )  * D1 

750  K=2*PI*F/C 

760  BH  =  A(2,I)/A  (1,1)  3  IH=A  (3 ,1) /A  (1 , 1) 

770  AA=SIN  (K*L)  3  B3=SI  N  (K*  (L-S)  )  3  CC=COS  (K*L)  2  DD=C0S 

(K*  (L-S)  ) 

780  IN  1 =  BB- AA*R  H  3  RN  1 =  AA  *IH 

790  RD1=RH*CC-DD  3  ID1=CC*IH 

800  M2  =  RD1  *RD1  +  ID1*ID1 

810  R(I)  ,RZ=  (RN1*RD1+IN  1*ID1)/M2  3  M  (I)  ,IZ=  (IN  1*RD  1-RN1* 

ID  1)  /«2 

820  RD=(RZ+1)/2 

330  ID=IZ/2  3  M=RD*RD+ID* TD 
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840  R?HO  =  -  (RD/M-1)  3  IP.H3  =~  ( ID/M) 

850  A 1  (I)  =SQR  (  RRHG*RRKO*I  SHO*IRHO)  3  PHI  ( I)  =-  AT  N  2 

87  0  next  I 

830  !  plot  results 

890  plotter  is  705  3  lorg  5  3  F0=5Q0  3  F1=1600 
900  Y 1  =-  1 0  3  Y 2=  1  0  3  S  1  =  1 

910  disp  "resistance"  3  1=1  3  L3="  RESISTIVE"  3  B3 

IMPEDANCE" 

920  gcsub  1070  3  if  Z1  then  1020 
930  Y 1 =- 1 0  3  Y 2= 1 0  3  S1=2 

940  disp  "reactance  3  1=2  3  L$="  REACTIVE"  3  B5=" 

IMPEDANCE" 

950  gosub  1070  3  if  Z1  then  1020 
960  Y 1  =0  3  Y2-1.4  3  SI  =.  2 

970  disp  "reflection"  3  1=3  3  LS="  REFLECTION"  3  E$= 

COEFFICIENT" 

980  gosub  1070  3  if  Z1  then  1020 
990  Y1=- 1 30  3  Y2=1  80  3  Sl=45 

1000  disp  "phase  angle"  3  1=4  3  L$="  PHASE"  3  £$=" 

ANGLE  (DEG)" 

1010  gcsub  1070  3  if  not  Z1  then  goto  430 
1020  disp  "lowest  Y"  ;3  input  Y1 
1030  disp  "highest  Y"  ;3  input  Y2 
1040  disp  "increments  of"  ;3  input  SI 
1050  or.  I  goto  91  0,940,970,1000 
1060  I  plotting  subroutines 

1070  disp  "change  paper  and  press  end  line" ; 3input  k$ 

1080  Y0=  (Y2-Y1) /7/2  5  3  Y3=  (Y2-Y1) /7*  1.  9  3  Y4=  ( Y 2- Y  1  0/7*1  .6 

3  Y5=Y1- Y 3  3  Y6=Y2+Y4 

1090  lorg  5  3  scale  -  (FI/5 )  ,  1 .  16  *F1 ,  Y5  ,  Y6  3  xaxis  Y1,FO,0, 
FI  3  yaxis  0,S1,Y1,Y2 
1100  Y7=Y  1 -Y3/6  3  Y8=Y1-Y3/2  3  ldir  0 
1110  for  X =0  to  FI  step  F0  3  move  X,Y7  3  label  VALS 
(X/1000)  3  next  X 

1120  move  .4*F1,Y8  3  label  "FREQUENCY-KHZ" 


1130 

1  la  i 

1150 
1160 
1170 
1180 
1190 
1200 
1210 
1220 

1230 
1240 
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for  Y=Y1  to  Y2  step  SI  >5  nove  -  (FI/20),  Y  3  xaoel  V”  \l~ 

(Y)  3  next  Y 

12  ir  PI/ 2  3  Y9  =  (Y2-Y1 ) /3  +  Y1  3  move  -(71/10), Y9  D  larel  • 

Li,BS  3  F9  =0 

for  L= 1  to  256  dtep  3 

on  I  goto  1170,1180,1190,1200 

F2=P.  (L)  3  goto  1210  • 

F2=M(L)  3  goto  1210 
F2=A1 (L)  3  goto  121  0 
F2=PHI (L) /PI  *  180 

move  F9,F2  3  F9  =  F9  +  3*D1  3  imove  -  (FI/500)  ,  Y0/2  • 

idraw  FI/250, 0  3  idraw  0,-Y0  3  idraw  -(FI/250),)  3 
idraw  0,Y0 
next  L 

dis?  ”1  for  redraw,  0  to  continue" ;3  input  zi  3  return  * 
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